Experimental scaling law for the sub-critical transition to turbulence in plane 

Poiseuille flow 



Gregoire Lemoultj^] Jean-Luc Aider and Jose Eduardo WesfreicQ 
PMMH (UMR7636 CNRS - ESPCI - UPMC Univ. Paris 6 - UPD Umv. Pans 7), 
Ecole Superieure de Physique et de Chimie Industnelles, 10 rue Vauquelin, 75005 Paris, France 

(Dated: December 1, 2011) 

We present an experimental study of transition to turbulence in a plane Poiseuille flow. Using a 
well-controlled perturbation, we analyse the flow using extensive Particule Image Velocimetry and 
flow visualisation (using Laser Induced Fluorescence) measurements and use the deformation of the 
mean velocity profile as a criterion to characterize the state of the flow. From a large parametric 
study, four different states are defined depending on the values of the Reynolds number and the 
amplitude of the perturbation. We discuss the role of coherent structures, like hairpin vortices, in 
the transition. We find that the minimal amplitude of the perturbation triggering transition scales 
like Re~^. 



PACS numbers: 47.27.Cn 

For more than a century [Ij, the transition to turbu- 
lence in shear flows has been a prohfic domain of study. 
Despite many theoretical investigations [2], it has not 
been possible to predict correctly this transitional pro- 
cess even for flows in simple geometries such as circular or 
plane Poiseuille flow, or plane Couette flow. Nowadays, 
the transition process in pipe and channel flows remains 
one of the most fundamental and practical problems still 
unsolved in fluid dynamics. 

Linear stability theory has been applied to plane 
Poiseuille flow, linearizing the Navier-Stokes equation 
near the stable parabolic proflle. The smallest unstable 
or critical Reynolds number obtained was Rcc = 5772 
(Re = Ucih/v with u^i the laminar center line velocity, 
h the half channel height and v the kinematic viscos- 
ity of the fluid) [3J |3]. This result is in contradiction 
with the experimental work of Carlson et al. who found 
1000 < Rec < 2000 0. 

Other approaches, related to the non-normal charac- 
ter of the Navier Stokes operator linearized around the 
stable laminar flow solution, can support important tran- 
sient growth of finite amplitude disturbances, related to 
streamwise and quasi-streamwise alignment [6l [7] . 

Waleffe, from an 3D non-linear modal reduction of the 
Navier-Stokes equations, adopted the idea of a self sus- 
tained process as origin of the transition [5] . After desta- 
bilization of these streaks, streamwise modes appear and 
regenerate the vortices through a non-linear interaction, 
as experimentally observed by Duriez et al. 9 . 

Transition to turbulence of sheared flows, is character- 
ized by a double threshold, in the initial disturbance and 
in the Reynolds number, where the larger the Reynolds 
number, the smaller the necessary perturbation. This be- 
havior, when the flow becomes unstable with respect to 
flnite amplitude disturbances, is described with a power 
law for the minimal amplitude of the disturbance trig- 
gering the transition: 

e = 0{Re^) 



After the first studies of Trefethen ^ on the critical 
exponent 7, Chapman studied the transient growth of 
sub-critical transition process in the plane Poiseuille flow 
[7]. He found 7 = —3/2 for initial streamwise vortices 
and 7 = —5/4 for initial oblique vortices. 

On the other hand, the results of the nonlinear study of 
Waleffe and Wang show that 7 = — 1 in shear flows |10| . 
Numerical experiments of Kreiss et al. give —21/4 < 7 < 

-1 m- 

Critical exponents were measured experimentally in 
plane Couette flow by Dauchot and Daviaud [12 , in pipes 
by Mullin et al. [13l [M] and in channel flow by Ben- 
Dov et al. [TS]. The latter found, using a cross jet as 
disturbance, that the critical jet velocity triggering the 
appearance of hairpin vortices scales like Re~^^^ . 

The present work studies experimentally the sub- 
critical transition of a plane Poiseuille flow perturbed by 
streamwise vortices induced by jets in cross flow. We will 
study the transition of the Poiseuille flow using the defor- 
mation of the mean velocity profile as a criterion instead 
of focusing on the apparition of hairpin vortices which is 
a step in the transition process. 

The mean flow distortion was described in channels 
by Eliahou et al. [TB]. It is indeed now well established, 
that the mean flow distortion is related to the presence of 
streamwise and quasi-streamwise elongated structures in 
the wall regions of boundary layers [17] , pipes and chan- 
nels. These structures, alternating low and high momen- 
tum contributions to the flow, generate through a nonlin- 
ear coupling, a global modification of the flow. Recently, 
Barkley gave a nearly complete model of the transition in 
pipe flow, introducing the modification of the mean ve- 
locity profile as one of the main ingredients of a system 
of coupled non-linear equations [T5] . 

When comparing experimental data to models, it is 
difficult to define properly the amplitude of the pertur- 
bation jT9] and also its critical value triggering the tran- 
sition. Our choice of selecting the mean fiow distortion 
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FIG. 1. Schematic view of the dedicated water channel used in PMMH. The developement section is 1 m long and the test 
section is 1.2 m long. Its cross section is 20 x 150 mm^. Left-bottom, scheme of the perturbation generation by continuous 
injection of water through small holes in the upper plate. Arrows represents the lift-up effect generating high and low speed 
streaks. Top-right, scheme of the multi plane PIV technique used to compute the average velocity profile. We used 11 planes 
separated by 3 mm. 



to characterize the transition puts a special emphasis on 
a more rigorous definition of the onset of turbulence. 

The experimental system is composed of a 3 m long 
plexiglass channel (Figure [I]). The test section's half 
height is h — 10mm, its length 220ft- and its width 15h. 
The perturbation is generated lOO/i downstream from the 
inlet to ensure a fully developed Poiseuille flow for all Re. 
The x,y and z axis are respectively the streamwise, nor- 
mal to the walls and spanwise coordinates, with y = 
in the middle of the channel and a: = where pertur- 
bations are injected. The design of the inlet section, to- 
gether with the smooth connections between all parts of 
the channel, minimize the upstream perturbations lead- 
ing to a laminar base flow until at least Re — 5500 (The 
Reynolds number is estimated from volume flow mea- 
surements) . 

The flow is perturbed by continuous injection of water 
through four circular holes, normal to the flow, drilled 
into the the upper wall with diameter d = 0.2h and spac- 
ing A = 3/i. The structure of the flow induced by the jets 
may be complex and depends strongly on the amplitude 
of the perturbation [201 HI] defined as the velocity ratio 
A = Ujet/uci, where Ujet is the mean jet velocity and Ud 
the unperturbed center line velocity . For < A < 2, 
one can consider that each jet creates a pair of counter- 
rotating streamwise vortices (figure [T]) similar to those 
created by solid vortex generators in a flat-plate bound- 



ary layer f9l. It has been shown that the destabilization 
of the streaks is a key step in a self-sustaining process be- 
tween streaks and streamwise vortices in Poiseuille flow 
as well as in a flat-plate boundary layer [21 [H]. It is 
also a key step in the transition scenario proposed by 
Chapman [7]. 
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FIG. 2. LIF visualization of the transition in the y = plane. 
The fiow goes from left to right. 

We show on figure [2] an exemple of visualization, ob- 
tained by Laser Induced Fluorescence (LIF) , of the tran- 
sition induced by the jets relatively far downstream from 
the injection. For x < 60h, one can see that the stream- 
wise vortices induced by the four jets are already unstable 
with clear streamwise modulation characteristic of hair- 
pin vortices. This correspond to the mixed behaviour 
observed by Tasaka et al. |23j . The transition to tur- 
bulence occurs further downstream [x > 60h) over the 
entire channel width. In the following, all measurements 
will be carried on in the region 75h < x < 80h illustrated 
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FIG. 3. From left to right, decomposition of the steps to compute the mean velocity profile. First we measure the time averaged 
(30 snapshots at AHz) velocity field in 11 a:: — j/ planes (only 6 presented here), then we average in the x-direction and finally, 
we average in the z-direction. 



on figure [T] 

The velocity field is studied using Particle Image Ve- 
locimetry (PIV) measurements ( figure [T] and [s]). The fluid 
is seeded with neutrally buoyant particles (dp « 5/im). 
To take into account the spanwise modulation of the 
flow, we measure velocity fields in a volume (Aa; = 3h, 
Az = A = 3/i, Ay = 2h), centred around x = 78h. The 
measurement volume is divided in llx — z planes between 
z = ±A/2 around a jet. For each plane, 30 instantaneous 
snapshots are taken at 4Hz giving 11 time-averaged ve- 
locity fields (figure [3ja). Then the PIV fields are space 
averaged over the streamwise direction (figure [3jb) . Fi- 
nally, the 11 profiles are averaged along the spanwise di- 
rection giving the final mean velocity profile (figure |3]c). 
Each mean profile presented in the following is the re- 
sult of an average over 26000 profiles. The total time to 
measure one mean profile is 3 minutes, thus averaging 
intermittent effects. 

The transition from the laminar to the turbulent 
regime is associated with a deformation of the mean ve- 
locity profile from a parabolic to a plug profile (figure |4]). 
The profile can then be used as a quantitative criterion 
to define the different states of the flow. We construct a 
state parameter u as: 

max(u) 
u = 

Ucl 

where u is the perturbed mean velocity profile. 

We plot on figure [5] the evolution of w as a function of 
Re (500 < Re < 5300) and A (0 < ^ < 2). The state 
plane is sampled with 450 points. The sampling has been 
refined in the transition region. Two different domains or 
plateaux can be clearly identified: the laminar (red) and 
turbulent (blue) state, separated by a sharp cliff, giving 
a description of the global state of the flow. We now need 
to define a rigorous criterion to identify the transition. 

For this purpose, we visualized the flow using LIF 
in the x = 80h plane. The structure of the flow is 
clearly modified when the amplitude of perturbation is 
increased. It is illustrated on figure |4] for Re = 2830. For 




y/h 

FIG. 4. First column represent the evolution of the mean ve- 
locity profiles, measured by PIV at a; = SOh, from which the 
parameter u is defined and in following columns two snap- 
shots taken using the LIF technique, at different times, in a 
y — z plane located at a:: = SOh. From top to bottom, rows 
correspond respectively to A = 0.276, 0.248, 0.209 and 0.167, 
Re = 2830 for all rows. 



a small perturbation (fourth row, A — 0.167), the flow 
is steady and laminar. Even if the fluorescent patches 
show traces of hairpin vortices, the mean velocity pro- 
file remains parabolic {u — I). As the amplitude of 
the perturbation is increased (third row, A = 0.209), 
hairpin vortices becomes larger but still smaller than 
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FIG. 5. 3D plot of the surface representing the parameter u = 
max(tt) /uci as a function of Re and A showing well separated 
domains of each state: laminar and turbulent. 



the channel's half-height and slightly unsteady. In this 
case the mean velocity profile becomes asymmetric and 
u decreases {u — 0.9). For A — 0.248 (second row), 
hairpin vortices becomes unstable, with intermittent ex- 
cursions in the lower half of the channel. The velocity 
profile becomes flatter, more symmetric and u still de- 
creases (m = 0.82). Finally, for A — 0.276 (first row), the 
flow becomes turbulent with large mixing over the entire 
cross section, leading to a typical turbulent plug proflle 
(m = 0.8). 

Thanks to this analysis, it is now possible to define a 
quantitative criterion for the onset of turbulence. In the 
following, we will consider that the flow becomes turbu- 
lent when u 0.8. For 0.8 < m < 1, the flow is in an 
intermediate state dominated by hairpin vortices. 



We plot on flgure [6] the state diagram of the flow, 
u{Re,A), as a contour plot in a log-log representation. 
We define four different regions representing different 
states of the flow. In region I, u = 1, the flow is steady 
and the mean velocity profile is parabolic. This region 
corresponds to the laminar state. Region II (u « 0.8) 
corresponds to the turbulent state, with high mixing and 
a plug velocity profile. Region III (0.8 ^ m < 1) cor- 
responds to a narrow intermediate state dominated by 
stable hairpin vortices, similar to those observed in jets 
in cross flow in a laminar boundary layer |21| . In region 
IV, perturbations induced by continuous jets do not trig- 
ger the transition but we observe some long-lived flows 
different from the well known parabolic proflle. 




FIG. 7. Log-log plot of the stability curve for the onset of 
turbulence in plane Poiseuille flow. Squares correspond to 
u = 0.81 and solid line is the -1 slope. The minimal ampli- 
tude triggering the transition scales like Re~^ at higher Re 
numbers. 
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FIG. 6. Contour plot in log-log scale of u{Re, A) showing 
the power law for the onset of turbulence for Re > 1500. I: 
laminar state, II: turbulent state. III: Regime dominated by 
hairpin vortices and IV: Below Re — 1500 it is not possible 
to sustain turbulence by this mean. 



For Re > 1500, we define the transition line as the 
separation between regions II and III, corresponding to 
u = 0.81. We show, in figure [7] a log-log plot of the min- 
imal amplitude of the perturbation Ac = {ujet/uci)c = 
f{Re) triggering the transition, defined as u{Ac) = 0.81. 
The —1 slope proposed by Waleffe is added on the plot, 
showing good agreement for Re > 2000. For Re < 2000, 
the asymptotic regime is not reached and experimental 
points deviate from the —1 slope. Using Ac for Re < 2000 
to evaluate 7, would lead to an underestimate of the ex- 
ponent [15] . 

The sub-critical transition of plane Poiseuille flow has 
been studied quantitatively through extensive PIV and 
LIF measurements. A well-defined state variable, u, has 
been introduced. Thanks to a large number of u measure- 
ments, we could define four different regimes of the flow. 
We then focused on the minimal amplitude of the per- 
turbation, Ac, triggering the transition. We found that 
Ac{Re) scales like Re~^ in agreement with the asymp- 
totic non linear theoretical model proposed by Waleffe 
for shear flows [TU]. We also show for the first time the 
role played by hairpin-like vortices in each step of the 
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transition. 
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